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Band structures and optical matrix elements of strained multiple quantum-wires (QWR's) are 
investigated theoretically via the effective bond-orbital model, which takes into account the effects 
of valence-band anisotropy and the band mixing. In particular, the Gai_2,In3;As QWR's grown 
by strain-induced lateral ordering (SILO) are considered. Recently, long wavelength Gai^ajIn^As 
QWR lasers have been fabricated via a single step molecular beam epitaxy technique which uses the 
SILO process. [1] Low threshold current and high optical anisotropy have been achieved. Multi-axial 
strains [combinations of biaxial strains in the (001) and (110) planes] for QWR's are considered. Our 
calculated anisotropy in optical matrix elements (for light polarized parallel versus perpendicular 
to the QWR's axis) is in good agreement with experiment. We also find that the strain tends to 
increase the quantum confinement and enhance the anisotropy of the optical transitions. 



I. INTRODUCTION 



GuxIrii-xAs is one of the most important ternary III-V compound semiconductors. Its bandgap covers both the 
1.3 and 1.55 /xm range, which are the preferred wavelengths in long distance fiber communications [1,2]. However, 
long-wavclcngth photonic devices based on lattice-matched Gao.47/no.53As/InP hctcrostriicturcs suffer from strong 
Auger recombination and intervalence band absorption proccsses[7-10]. Recently, to improve the performance of long- 
wavelength semiconductor lasers, long wavelength('~ 1.55/xm) Ga^Ini-xAs quantum-wire (QWR) lasers have been 
grown by a single step molecular beam epitaxy techniquc[l]. It is found that the QWR laser structures are a promising 
choice because of the many predicted benefits, such as higher gain, reduced temperature sensitivity, higher modulation 
bandwidths, and narrower spectral linewidths[2]. Short-period superlattices (SPS) in the direction perpendicular to 
the growth direction can be formed via the strain-induced lateral-layer ordering (SILO) process[l,2]. A quantum wire 
heterostructure can then be created by simply utilizing this SPS structure as the quantum well region in a conventional 
quantum heterostructure. Besides the self- assembled lateral ordering, it is believed that the strain also plays a key 
role [2] in the temperature stability and high anisotropy for the QWR laser structure. Recent studies indicate that 
the use of the strained-layer quantum wire heterostructure has advantages of high-quality interfaces and band-gap 
tuning which are independent of the lattice constant of the constituent materials[3,4,6]. Furthermore, much work has 
been undertaken which predicts that by using strained-layer superlattice to form the active region of a quantum-wire 
laser, the threshold current can be decreased by one order of magnitude, and the optical loss due to intervalence- 
band absorption and Auger recombination will also be greatly reduced[l,2,7-10]. Also, the temperature sensitivity is 
reduced by a order of magnitude compared with strain-free structures [2]. A typical temperature sensitivity of the lasing 
wavelength is(~ 5A/°C) for the usual GalnAs/lnP lasers. By using a distributed-feedback structure, the temperature 
dependence of the lasing wavelength is reduced to lA/°C[2]. With this strained GalnAs QWR, the dependence is 
smaller than 0.lA/''C[2]. It should be an important improvement on the current technology in fabricating the long 
wavelength lasers for fiber communication. In this paper we study the effects of multi-axial strain on the electronic 
and optical properties of the QWR structures grown via the SILO process. Through this study, we can gain a better 
understanding of the strain engineering of QWR structures suitable for the application of fiber-optical communication. 

In the experiment performed by Chou et al[l], the QWR active region is created in situ by the SILO process within 
the {GaAs)2/ {InAs)2.2 SPS regions. The SILO process generates a strong Ga/In lateral composition modulation 
and creates In-rich GaxIn\-xAs lateral QWs in the [110] direction. By sandwiching the In composition modulated 
layer (with x varying from 0.3 to 0.7) between A^o.24G'ao.24J"?^o.52As barrier layers, a strained QWR heterostructure 
is formcd[l,2]. The model QWR structure considered in the present paper is depicted in Figure. 1. The central region 
consists of the In composition modulated GaxIni-xAs layer (A/B/A/B...), in which the A strip is Ga-rich while 
the B strip is In-rich. The lateral period is L2 = L2 + if. In our calculation, we consider = Xf = 50 A and 
if = LY = 50A. We consider two cases of composition modulation. In case one, the Ga composition changes abruptly 
with X ~ 0.6(or0.7) in strip A and x = 0.4(or0.3) in strip B. In case two, x varies sinusoidally from 0.6 (or 0.7) at the 
center of strip A to 0.4 (or 0.3) at the center of strip B (with x = 0.5 at the border between A and B strips). 

II. THEORETICAL APPROACH 

The method used in this paper for calculating the strained quantum well band structure is based on the effective 
bond-orbital model (EBOM). A detailed description of this method without the effect of strain has been published 
elsewhere [5]. Here we give a brief account of the method and discuss the effect of strains. The effective bond-orbital 
model employed here is a tight-binding-like model which includes nearest-neighbor interactions among bond orbital 
residing on an fee lattice. The version of the EBOM used in this calculations describes the coupling between the upper 
four spin-3/2 valence bands and the lowest two s-like conduction bands. For our QWR system, the split-off bands 
are ignored since the spin-orbital splitting is large. A bond orbital is defined to be the proper linear combination of 
the two atomic orbitals within a zinc-blende crystal which best describes the states near the center of the Brillouin 
zone. The parameters that appear in the theory are given by a correspondence which is made by requiring that the 
Hamiltonian in the bond-orbital basis, when written for the bulk material and expanded to the second order in k, 
agree with the Luttinger-Kohn expression. For interactions across heterojunctions, we take the average of the matrix 
elements for the two bulk materials. EBOM has been successfully applied to the calculation of electronic states of 
superlattices [5], strained-layer superlattices [6,19], semiconductor quantum wires[15] and dots[17], and impurity levels 
in quantum dots[17] and quantum wells[17]. Here, we apply it to this multi-axially strained QWR structures. 

We solve the QWR electronic states according to a procedure similar to that described in [15]. We first calcu- 
late the band structure of a superlattice(SL) with M bilayers of GaxIni-xAs (the well material) and N bilayers of 
Alo.24Gao,24lno,52As (the barrier material) using the slab method. The superlattice period is L± = + L^= 



(M+N)a/2, where a is the lattice constant. For a sufficiently large value of the barrier thickness (we used = 50 A) 
the QW's arc essentially decoupled and wc ignore the dependence on the wave number in the growth direction, taking 
kz = 0. The SL eigenstates with fc^ = arc denoted by \ki,k2,n >, i.e, 

H^^\ki,k2,n>= En{ki,k2)\ki,k2,n> (1) 

where H^^ is the Hamiltonian for the Guylni-y As / Alf).24Gao.24lno. 52^3 SL (y = 0.6or0.7) , where fci and k2 are the 
two components of the SL in-plane wave vector along and perpendicular to the QWR axis, respectively, fci is a good 
quantum number for the QWR array, while ^2 is not. The QWR wave function is a linear combination of SL wave 
functions with the same ki but with values of ^2 differing by 2tt/L2. We define k2 to be the wave number of the 
QWR array in the [110] direction, and we can write ^2 = ^2 + 27rC/i2 with ^ = 0, 1, 2, ...,m, where m = £2/02 with 
a2 = a/V2. The QWR Hamihonian is ijQ^-R = H^'^ + H', where H' = Hoa^im-^As - Hoayim-yAs is treated as a 
perturbation. The Hamiltonian matrix for QWR 

M(n, C|n', CO =< kuk2 + 2ttC,' /L2\H'^^'^\kx, fcs + 27rC/L2, n > (2) 

is then diagonalized in the basis consisting of SL eigenstates. The basis is truncated to reduce the computer time, 
thus making the method a variational one. The subbands closest to the band edge are converged to within 0.1 mev. 
The size of the Hamiltonian matrix, originally given by the number of sites in a QWR supercell times the number 
of coupled bands, is now reduced to the number of SL states necessary to give accurate results for the QWR of 
concern. In our case, this number needed to get well converged results for the two uppermost pair of valence subbands 
is around 200. However, we have used ~ 800 basis states in order to ensure the accuracy of the deeper lying sates 
for all calculations in this present paper. The effect of lateral intermixing of Ga and In is accounted for by using 
the virtual-crystal approximation. Bond-orbital parameters for the inter-difFuscd materials are obtained by linearly 
interpolating between the values of the parameters for the pure materials. The optical parameters in these expressions 
are determined by requiring that the optical matrix elements between bulk states obtained by EBOM be identical to 
those obtained in the k • p theory up to second order in k[15]. 

The effect of strain is included by adding a strain Hamiltonian iJ** to the EBOM effective Hamiltonian[6,10,16]. 
The matrix elements of H^* in the bond-orbital basis can be obtained by the deformation-potential theory of Bir 
and Pikus[16]. Here, we consider the combination of two bi-axial strains, one in the (001) plane (due to the lattice 
mismatch at the interfaces between GalnAs and AlGalnAs) and the other in the (110) plane (due to the Ga composition 
modulation ). The resulting strains in the SL layers consist of both hydrostatic and uniaxial components. For case 1 
(Ga composition changes abruptly from region A to region B), the lattice constant of the strained layers in the [001] 
direction (aooi) is given by minimizing the free energy of the fee system due to strain[ll]: 

P = 2^ll{^\lA + '^22A + <^33a)^^ (3) 
+ C^2{^11A^22A + ^22A^diA + ^11A^33a)L2 
+ 2^n(^llB + ^^223 + *^33B)-^f 
+ Ci2{illB^22B + ^228^338 + ^11B^33b)L2 , 

where (-^f ) is the lateral layer thickness and eijAiUjs) is the strain tensor in the Ga-rich (In-rich) region. C'^'s 
(C^'s) are the elastic constants for the Ga-rich (In-rich) materials. Here we have used a rotated Cartesian coordinates 
in which x' = [ll0],y' = [110], and z' — [001]. A constraint a"^ — = ajnP = cLAUnGaAs has been imposed to 
keep the A/q. 24^09. 24-/^ '^0.52^* barrier strain free, since it is lattice-matched to the InP substratc[l,2], where ajj* is the 
in-plane strained lattice constant (perpendicular to [001] axis) for material a. The multi-axial strain caused by the 
lattice-mismatch in both (001) and (110) planes is simply 

eila = e22a = , (4) 

CLct 

where oa (ob) is the unstrained lattice constant of GaxIni-xAs {Gai-xIUxAs) [x = 0.6 for case 1) and 

aOOl — da 

essa = (5) 



with efj = for i ^ j. 



The minimization procedure leads to 



OOOI = — : 



I'^A^ ^11^2 I '^B 



(6) 



which is the strained lattice constant in the [001] direction. 

Due to symmetry, the multi-axial strain tensor is diagonal in the x'y'z' coordinates, which leads to a diagonal strain 
Hamiltonian 





-IWh + D2 

-AVh + Ds 



(7) 



where 

AVh = (ai + a2)(eii + £22 + £33), -Di = &(2eii - £22 - £33), -D2 = &(2e22 - en - £33), -03 = 6(2e33 - £22 - en). 
The strain potential on the s states is given by [6] 

AVc = ci(en + £22 + ess), 

The strain Hamiltonian in the bond-orbital basis \JM > can be easily found by using the coupling constants [5], i.e, 

< JM\H'^\J'M' >= J2 C{a,a;J,MyCia',a;J',M')H'J^, (8) 



The clastic constants C12 and Cn for GaAs, InAs and AlAs can be found in Rof. [12,18]. The deformation 
potentialsai, 02, b, Ci can be found in Ref.[20,21]. The linear interpolation and virtual crystal approximation is 
used to obtain the corresponding parameters for the GalnAs and AlGalnAs. 

After rotating back to the original Cartesian coordinates, the final results for the strain Hamiltonian in the bond- 
orbital basis \JM > is 
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(9) 



with 



A = 



£)3 + 0.5(Z)i + ^2) 



B 



P1+D2) ^ ^ 2.5{Di+D2) 

L>3, O = 



One can see non-diagonal terms arise due to the rotation, which will cause additional mixing and optic transitions 
between different valence bands. 

The above strain Hamiltonian is derived for case one only, where the Ga composition (x) changes abruptly from 
region A to region B. For case two, x varies continuously from Xm to 1 —Xm {x„i = 0.6or0.7). In this case, wc shall first 
calculate the strain Hamiltonian for the abrupt case to get i?*^ and Hf_^^, using Eq. (8). We then obtain the strain 
Hamiltonian at any atomic layer with Ga composition x via a linear interpolation between and Hf_^^ . This is 
consistent in spirit with the virtual crystal approximation we have used for obtaining the interaction parameters in 
EBOM. 



III. RESULTS AND DISCUSSIONS 



In Figure 2, we show a schematic diagram indicating the alignment between band edges of the constituent materials 
for QWR in case 1 (square profile) with and without the effect of strain. All energies are measured with respect to the 
bulk InAs valence band edge. This diagram is useful in understanding the quantum confinement effect on the QWRs 



considered here. From Figure 2, we notice that the strain [in particular, the one in the (110) plane] has a significant 
effect on the band alignment. Take Figure 2(a) for example, without strain (solid lines) the conduction (valence) band 
of Gao.4Ino.6As is below (above) that of Gao.6l1io.4As by 247 meV (14 meV). With strain (dotted lines), the situation 
is reversed and the conduction (valence) band of Gao.4Ino.6As is above (below) that of Gao.6lno.4As by 312 meV (76 
meV); thus, both electrons and holes will be confined in the Ga-rich region of the QWR. We found that about 80 
% of these changes are caused by the strain in the (110) plane, since there is a stronger lattice mismatch between 
Gao.6lno.4As and Gao.4Ino.6As compared with that between Gao.6lno.4As and Alo.24Gao.24Ino.52As (or InP). 

The strain also causes a splitting between the heavy-hole (HH) and light-hole (LH) bands with the LH band lying 
above the HH band. However, the splitting is rather small, about 3 meV in Gao.6lno.4As. When the confinement 
effect due to the barrier material (Alo.24Gao.24Ino.52As) is included, the HH band again hes above the LH band, due 
the difference in the effective masses along the growth direction (z). The dash-dotted lines in Figure 2(a) indicate the 
superlattice band edges of 50AGao.6lno.4As (or Gao.4Ino.6As) sandwiched between 50AAlo.24Gao.24Ino.52As barriers. 
The difference in superlattice band edges (dash-dotted lines) between Gao.6lno.4As and Gao.4Ino.6As determines the 
degree of lateral quantum confinement in the QWR. For the present case, the conduction-band offset is 252 meV and 
the valence-band (for HH only) offset is 52 meV, as far as the lateral confinement is concerned. Both offsets are large 
enough to give rise to strong lateral confinement for electrons and holes in the Ga-rich region. 

Figure 3 shows the near zone-center valence subband structures of square-shaped QWRs with the multi-axial strain 
with (a) Ga composition (x) changing abruptly from 0.6 in region A to 0.4 in region B and (b)Ga composition (x) 
changing abruptly from 0.7 in region A to 0.3 in region B. To compare with experiment [1,2], we choose material param- 
eters appropriate for temperature at 77K[12]. All subbands are two-fold degenerate due to the Kramer's degeneracy 
and they are labeled according to the characters of their underlying Bloch functions at the zone-center: HH for heavy 
hole and LH for light-hole (only the confined subbands are labeled). The conduction subbands are approximately 
parabolic as usual with a zone-center subband minimum equal to 826 meV in (a) and 707 meV in (b) (not shown in 
the figure). This gives an energy gap 860 meV for the (x=0.6/0.4) QWR array and 708 meV for the (x=0. 7/0.3) QWR 
array . 

The observed Cl-HHl excitonic transition is at 735 meV (or 1.65/im) for the QWR[1,2] at 77K. In comparing the 
band gaps with the experiment, one should also take into account the exciton binding energy which is around 20 
meV for this size of QWR. Thus, the theoretical result for the (x=0.7/0.3) QWR array is in closer agreement with 
experiment, but about 50 meV too low. 

Comparing the band structures in both ki ([110]) and ^2 ([110]) directions, we noticed an apparent anisotropy in the 
energy dispersion. The dispersions in the k2 direction for the five uppermost (confined) valence bands are rather small, 
indicating strong lateral confinement. We observe strong anti-CjCOSsing effect between the HH1/HH2 and LHl subband 
at ki near 0.02A~^ similar to what happens in a quantum wellEj. Here the HH2 subband of the QWR corresponds to 
the [110] zone-folded part of the HHl subband of the quantum well ( the envelope function has odd parity in the [110] 
direction but even parity in the [110] direction). 

The subbands in Figure 3(b) have less dispersion in the ^2 direction compared with those in Figure 3(a) as a result 
of stronger lateral confinement. This is caused by the larger band discontinuities in the x=0.3/0.7 case (versus the 
x=0.4/0.6 case) as can be seen by comparing band alignments shown in Figure 2. 

Figure 4 shows the near zone-center valence subband structures of a QWR with a sinusoidal lateral modulation with 
Ga composition (x) varying as the position (y') in the [110] direction between two extreme values Xm and 1 — Xm, i.e. 

x{y') = 0.5 -I- (0.5 - Xra)sin{2TTy'/L2) 

for Xrn = 0.4 and (b) Xm — 0.3. In this case, the Hamiltonian including the multi-axial strain is calculated with a 
linear extroplation between the Hamiltonians for the maximum and minimum Ga composition, i.e. 

^Ga^Ini-^As ^ f Hca^^Im^^^As + (1 ^ f)HGai-^^In^^As, 

where Xm — 0.4 or 0.3 and / is determined by comparing the Ga composition on both sides of the equation 

X = X,nf + (1 - X,n)il - /), 

01 f = {1-Xm- X)/{1 - 2Xm)- 

The conduction band minimum (not shown) is 850 meV (757 meV) for the case Xm = 0.4 {xm = 0.3), which 
corresponds to a band gap of 889 meV (767 meV). Subtracting the exciton binding energy (~ 20 meV) from the band 
gap, we found that the Xm = 0.3 case is in fairly good agreement with the observed excitonic transition energy of 
735 meV. The band gaps of QWR's with the sinusoidal profile are consistently larger (by about 30-50 meV) than 
the corresponding QWR's with the square-profile (with the same extreme values of Ga composition). The sinusoidal 



variation in Ga composition gives rise to more energy, since the well region contains more materials with band gap 

higher than the minimum in comparison with the square-profile. Furthermore, the lateral strain is reduced (less 
mismatch on the average) which also tend to increase the band gap. 

Comparing Figure 4 with Figure 3, we notice that the spacing between HHl and HH2 is substantially larger for 
the sinusoidal profile than for the square profile. This can be understood by the following argument. The envelope 
function for the HH2 state is more spread out than that for the HHl state, thus its energy is increased more in the 
sinusoidal profile (with higher probability being in materials with higher band gap) compared with the square profile. 

Next we discuss the optical properties of QWR's. Since the QWR states can be qualitatively viewed as the zone- 
folded states of the superlattice (SL) states (the zeroth-order states) via the lateral confinement, it is instructive 
to examine the optical matrix elements of the SL case. Figure 5 shows the squared optical matrix elements of the 
(50A/50A) Gao.7Ino.3As/Alo.24Gao.24Ino.52As superlattice for the HHl-Cl and LHl-Cl transitions versus the wave 
vector k2- The solid and dashed lines are for the polarization vector along the fci ([110]) and ([110]) directions, 
respectively. We note that the optical matrix elements are isotropic in the x-y plane at the zone-center, while they 
become anisotropic at finite values of fe2- At the zone center, the HH state consists of Bloch states with atomic 
character (x' -I- iy') t (l-^,-^ >= |3/2, 3/2 >), while the LH state consists of Bloch states with atomic character 
(x' — iy') t +2z I (|J, M >= |3/2,— 1/2 >); thus, the corresponding optical transitions (to an s-like conduction band 
state) arc isotropic in the x-y plane. Here x'{y') is the coordinate along the fci(fc2) direction. For finite fc2, the HH 
and LH characters are mixed, with the HH state consisting of more .x'-character than y'-character, thus in favor of the 
polarization vector parallel to the fci direction. Note that the HH (LH) band tends to have an atomic character with 
polarization perpendicular (parallel) to the direction of wave vector. This is a direct consequence of the fact that the 
(pPt:) interaction is weaker that the (ppc) interaction in a tight-binding model. Since QWR states are derived from 
the SL states with finite values of ^2, we expect the optical matrix elements of QWR to be anisotropic in the x-y plane 
as well. 

Figures 6 and 7 show the squared optical matrix elements versus wave vectors of the corresponding QWR's considered 
in Figures 3 and 4, respectively, respectively. The solid and dashed lines are for the polarization vector along the 
ki ([110]) (parallel to the wire) and ^2 ([110]) (perpendicular to the wire) directions. The optical matrix elements 
together with the subband structures discussed above provide the essential ingredients for understanding the optical 
transitions observed in the photoluminescence (PL) measurements. In all figures the optical transitions considered are 
from the topmost three valence subbands to the lowest conduction subband: HHl-Cl, HH2-C1, and LHl-Cl. These 
curves were obtained by summing the contributions of two degenerate subbands (due to Kramer's degeneracy) for the 
initial and final states. To compare with experimental (PL) results, we concentrate on the HHl-Cl transition. 

From Figures 6 and 7, we found for the parallel polarization (solid lines) , the squared optical matrix clement for the 
HHl-Cl transition has a maximum at the zone center with a value near 25 eV and remains close to this maximum 
value for all finite ^2 and with fci = 0. For the perpendicular polarization, the HHl-Cl transition is very small for finite 
^2 and ki = 0. This means that the wave fimction of the HHl state has mostly x' character (and some z character). 
Thus, we conclude that the strong lateral confinement forces the wave function of the HHl state with fci = to change 
from the x' + iy' character (in the SL case) into mostly x' character. At finite fci, the character of the HHl state 
gradually changes back toward the x' + iy' character and the optical transition becomes almost isotropic when ki is 
comparable to zone-boundary value of k^iir / L2), and finally turns into mostly y' character as ki becomes much larger 
than (77/^2)- The anti-crossing behavior of the HH1/HH2 subbands with the LHl subband further complicates the 
whole picture at kx « 2A~^. 

The HH2-C1 transition has zero optical strength at the zone center for both polarizations. This is expected, since 

the HH2 state has odd parity in the envelope function, which leads to forbidden transition at the zone center. The 
symmetry restriction is relaxed as the wave vector deviate from zero. The LHl-Cl transition has large optical strength 
(around 17 eV) for the perpendicular polarization and very weak strength (around 1 eV) for the parallel polarization at 
the zone center. This indicates that the LHl states consists of mostly y' character. So the strong lateral confinement 
forces the the wave function of the LHl state with fci = to change from the x' — iy' character (in the SL case) into 
mostly y' character. In other words, the HHl state (with pure x' -\- iy' character) and LHl state (with pure x' — iy' 
character) in the SL case are mixed thoroughly by the lateral confinement in QWR to produce a predominantly x' 
state (HHl) and a predominantly y' (LHl) state. Note that all thse states consist of appreciable z character, which 
will appear in the optical transition with z polarization. 

To calculate the anisotropy, we integrate the squared optical matrix element over range of fci corresponding to 
the spread of exciton envelope function in the fci space. The exciton envelope funtion is obtained by solving the ID 
Schroding equation for the exciton in the effective-mass approximation 

[-|^(^)' + yxix')]F{x') = ExF{x'), (10) 



where is the exciton reduced mass (w 0.037mo) and Vx{x') is the effective ID exciton potential given b; 



The parameter /3 is obtained by extropolation from the values given in Table V of Ref.S. For a 50A quantum wire, we 
obtain f3 — 0.07A~^. Eq. (10) is solved numerically, and the exciton envelope function in ki space is obtained via the 
Fourier transform of F{x'). The exciton binding energy obtained is 23 meV. We found that the ratio of the averaged 
optical strength for the HHl-Cl transition for the parallel to perpendicular component of the polarization vector is 
2.65 (4.04) for the x = 0.4/0.6 {x = 0.3/0.7) QWR with square profile, and 2.86 (7.19) for the a;,„ = 0.4 (x„ = 0.3) 
QWR with sinusoidal profile. The Xm = 0.3 case shows stronger optical anisotropy than the Xm = 0.4 case, indicating 
that the stronger lateral confinement leads to stronger optical anisotropy as expected. Experimentally, the optical 
anisotropy is found to be around 2-4. So, our results are consistent with experiment. 



IV. CONCLUSION 



We have calculated the band structures and optical matrix elements for the strained QWR grown by the SILO 
method. The effect of multi-axial strain on the valence subbands and optical matrix elements is discussed. Our 
theoretical studies provide the explanation of the anisotropy in optical matrix elements of these QWR's observed 
experimentally. We find that the biaxial strain due to the lattice mismatch between the Ga-rich and In-rich regions is 
most dominant. It tends to increase the lateral confinement in the QWR and enhances the anisotropy of the optical 
transitions which may be useful for certain applications in optical communication. We also calculated the effect of 
lateral composition modulation on the band structures and optical properties and find that it increases the band gaps 
and reduces the optical anisotropy. From the above discussions, we can use the theoretical predictions to guide the 
engineering design of QWR optical devices. The temperature effect will be incorporated in our future studies in order 
to understand the temperature stability of the optical transitions of QWR's grown by the SILO method. 
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Figure Captions 

Fig. 1. Schematic sketch of the QWR array fabricated in Refs. [1,2]. The QWR axis hes in the [110] direction. A and 
B label the Ga-rich and In-rich strips, respectively. 

Fig. 2. Schematic diagram indicating the alignment between band edges of the constiuent materials for QWR in case 
1 (square profile) with and without the effect of strain. Solid lines are for unstrained bulk, dotted lines are for bulk 
under multi-axial strain appropriate for the present QWR, and dash-dotted lines are for (50A/50A) GalnAs/AlGalnAs 
superlattice including the multi-axial strain. 

Fig. 3. Valance subband structures for QWR's with square profile for (a) Ga-composition x = 0.6 in Ga-rich and 
X = 0.4 in In-rich region and (b) x = 0.7 in Ga-rich and x = 0.3 in In-rich region. 

Fig. 4. Valance subband structures for QWR's with sinusodial profile with x ranging from Xm to 1 — Xm for (a) 

■Tm = 0.6 and (b) Xm = 0.7. 

Fig. 5. Squared optical matrix elements for transitions from HHl and LHl to the first conduction subband of the 
(50A/50A) Gao.7Ino.3As/Alo.24Gao.24Ino. 52 As superlattice for light polarized parallel (solid), perpendicular (dashed) 
to the QWR axis and z component (dotted). 

Fig. 6. Squared optical matrix elements for transitions from the top three valence subbands to the first conduction 
subband for light polarized parallel (solid) and perpendicular (dashed) to the QWR axis for QWR's considered in Fig. 
3. 

Fig. 7. Squared optical matrix elements for transitions from the top three valence subbands to the first conduction 
subband for light polarized parallel (solid) and perpendicular (dashed) to the QWR axis for QWR's considered in Fig. 
4. 
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